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We use an effective Markovian description to study the long-time behaviour of a nonlinear second 
order Langevin equation with Gaussian noise. When dissipation is neglected, the energy of the 
system grows as with time a power-law with an anomalous scaling exponent that depends both 
on the confining potential and on the high frequency distribution of the noise. The asymptotic 
expression of the Probability Distribution Function in phase space is calculated analytically. The 
results are extended to the case where small dissipative effects are taken into account. 
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The influence of a random perturbation on a dynamical system is a problem of interest in various fields of science and 
engineering [2, d, 0, 0, S @] • The first example of a differential equation with stochastic terms appeared in Langevin's 
CTS ' study of Brownian motion 0, H[ : Langevin modelized the action of the solvant molecules on the Brownian particle 
as the sum of a deterministic viscous friction, proportional to the velocity of the Brownian particle, and of a random 
force of auto-correlation proportional to the temperature of the bath. Since then, it has been customary to add in the 
dynamical equations some phenomenological stochastic terms that describe random environmental loadings (e.g. the 
influence of a turbulent wind on a suspension bridge, or the study of random parametric vibration of helicopter rotor 
blades in atmospheric turbulent flow [13] ). Of particular interest is the determination of the energy flow into the 
system from external sources when the characteristic time of the parameters variations matches one of the natural 
frequencies of the system: parametric resonance then occurs and the rate of increase of the amplitude is generally 
OA . exponential leading to an instability. The growth of the response is limited by various nonlinear effects. 

Several methods have been developed to study random parametric vibrations [T^. One of the most efficient 
techniques is the averaging principle developed by Bogoliubov and Mitropol'skii for deterministic nonlinear vibrations 
Til ] where rapidly fluctuating circular coordinates are averaged out leading to a set of effective dynamical equations 



< 



for slow variables. This method was extended to stochastic systems by Stratonovich [6(] and put on a rigorous 
mathematical basis by Khas'minskii [ijj and by Papanicolaou and Kohler [13| . Since then, stochastic averaging has 
become a powerful method 14l_ 15| ( for a recent review see [16| and references therein). 
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In a series of recent works [17|, |18|, |19|, l20(, we studied the long time behaviour of the nonlinear oscillator subject 
to parametric noise. We showed that, in the absence of dissipation, the nonlinear terms in the potential stiffness 
I ] inhibit the exponential growth of the amplitude. The observables of the system (the amplitude, the momentum, 
■ the energy) rather display power-law scalings with anomalous diffusion exponents. When the parametric noise is a 
t-H | Gaussian white noi se ( ».e. , it has a vanishingly small correlation-time), the averaging method, applied to the energy 
envelope @, [l(| HH HH, [H| , allows us to calculate analytically the time-asymptotic probability distribution function 
(PDF) of the syste m in phase space knowing this PDF, the scaling exponents and the corresponding prefactors 
are readily deduced. However, for colored noise a competition between conflicting time scales occurs. In fact, the 
nonlinear oscillator has an amplitude-dependent intrinsic frequency that increases with the amplitude. And as the 
oscillator absorbs energy from its environment, its amplitude grows and at a certain stage, the intrinsic period becomes 
smaller than the correlation time of the noise. This corresponds to a cross-over regime at which the correlation time 
of the noise ceases to be the smallest time-scale in the system: the scaling laws that govern the growth of amplitude, 
momentum and energy differ from those calculated for a white noise. Thus, when the amplitude of the oscillator is 
small its intrinsic period is large and the noise appears as if it were white and white noise exponents prevail; but at 
large amplitudes, the scaling regime changes and new exponents appear. Because of these conflicting time scales, the 
averaging technique is difficult to implement for a colored noise: at lowest order the noise itself is averaged out and 
the energy transfer stops at the cross-over time. Therefore one has to perform averaging at higher orders. When the 
colored noise is an Ornstein-Uhlenbeck process the calculations can be carried out by a second-order averaging, which 
requires rather tedious mathematical manipulations [l8L [l9l ] . It is also possible to calculate the cross-over between 
the white noise and the Ornstein-Uhlenbeck scaling regimes. The averaging method works at second order for the 
Ornstein-Uhlenbeck noise because its time derivative is a white noise. However, if the random excitation is generated 
from a white noise through a differential equation of order n, one has to perform averaging at (n + l)-th order and in 
practice, the calculations are intractable. 

In the present work, we follow an entirely different approach to study the nonlinear oscillator subject to a parametric 
Gaussian noise with an arbitrary spectrum (with the assumption that the spectrum decays as a power-law at high 
frequencies). We shall use an effective coarse-grained Markovian description of the dynamics, following a technique 
developed by Carmeli and Nitzan [24| (see also [23|,[25| for a similar approach). This technique will allow us to calculate 



analytically the asymptotic PDF which leads to the formulae for the growth of the amplitude, of the momentum and 
the energy transfer. In particular we shall prove that the scaling exponents depend both on the stiffness of the 
potential at infinity and on the smoothness of the random excitation; smoother the noise (which corresponds to a 
faster decay of the power spectrum at high frequencies), less efficient is the energy transfer from the bath to the 
oscillator. The method used here can be adapted both to additive and to multiplicative noise, and can also be used 
when a small friction is present: the system reaches at large times a non-equilibrium steady state in which physical 
observables do not grow anymore; the cross-over from power-law growth to this steady state occurs when the rate of 
energy dissipation by friction matches that of energy absorption from the random environmental loading. 

The outline of this work is as follows. In Section [H we define precisely the model we shall study. In Section [II] we 
use the underlying integrability of the system to write exact dynamical equations in energy-angle variables and we 
use the coarse-grained Markovian description to derive an effective Fokker-Planck equation for the energy variable. In 
Section Hill we derive explicit formulae for various cases: multiplicative or additive noise, with or without dissipation. 
This leads to a rather exhaustive description of all the different cases. In particular, we verify that this method allows 
us to recover the analytical results obtained previously for white and Ornstein-Uhlenbeck noises. The last section is 
devoted to concluding remarks. 



I. THE NONLINEAR OSCILLATOR WITH PARAMETRIC NOISE 



A paradigm for the study of interplay of noise and nonlinearity is the nonlinear oscillator subject to parametric 
random excitations: 

*L X (t) + 7 fU w + (a ,g + m X (t) + ^ = o . (i) 

The variable x(t) represents the amplitude of the oscillator at time t. The potential U(x) that confines the oscillator 
is assumed to grow faster than quadratically when \x\ — > oo giving rise to nonlinear terms in the restoring force. We 
shall make the simplifying assumption that IA is an even function of x and behaves as a power-law of x when |i| — > oo. 
Then, a suitable rescaling of x allows us to write 

U with v > 2 . (2) 

Typically, v is an integer; the value v = 3 corresponds to the Duffing oscillator. 

The physical interpretation of Eq. ([T]) is that the linear stiffness of the oscillator fluctuates around its mean value 
u>q because of randomness in the external conditions and this randomness is represented by external noise We 
also suppose that the oscillator is subject to a linear friction with damping coefficient 7. 

Equation (JTJ) is thus a nonlinear stochastic differential equation. When the multiplicative noise £(t) is a white 
noise, a coherent convention to perform stochastic calculus must be chosen. Although Ito calculus is favored by 
mathematicians, we shall use here the Stratonovich calculus [HI, Ht| which is physically more sound because it 
appears naturally when one considers the white noise as a limit of colored noise with very short correlation time 
0, dfl . This equation seems to be very elementary but it embodies many features of random dynamics: inertial 
effects, nonlinear stiffness and parametric noise. In fact, many complex dynamical systems that appear in realistic 
engineering problems can be reduced after some simplifying assumptions to Equation |T]). For example, the torsional 
stability of a suspension bridge under the influence of wind loads can be reduced to an equation similar to the one we 
are studying [13, HI] ; similarly, the dynamics of liquid sloshing, the roll motion of a ship or the stability of helicopter 
rotor blades in hoover flight under atmospheric turbulence can be reduced to effective single-degree-of-freedom systems 
represented by a second order equation with random parametric vibrations (for explicit derivation of such equations 
see e.g. fr(|). A similar equation has also been proposed by Fermi as an acceleration mechanism for interstellar 
particles [29l l30l . l3lj ] . Finally, from the mathematical point of view, Equation ([T]) is also very appealing: it is rich 
enough to exhibit an interesting dynamical behaviour but simple enough to allow for explicit solutions [32L I33L H3] . 
This explains why such a simple model can play the role of a paradigm. 

The phase space origin, x = and dx/dt = is a solution of Eq. (JTJ) . However, it can be shown that this solution is 
unstable [H, [36| when the power-spectrum of the noise contains all possible frequencies. When friction is neglected 
(i.e. if the underlying deterministic system is Hamiltonian) , then because of the permanent injection of energy into 
the system by the noise, the amplitude, the velocity and the energy undergo anomalous diffusion. The associated 
anomalous diffusion exponents and amplitudes have been calculated exactly when the random excitation is a Gaussian 
white noise [T3| or an Ornstein-Uhlenbeck process 

In the present work, we study the effect of the statistical properties of £(i) on the long time behaviour of the 
dynamical variable x(t). We must therefore specify the characteristics of the random perturbation We shall 



consider a stationary Gaussian noise of zero mean value. A Gaussian process is fully characterized by its auto- 
correlation function defined as 



Stf-t) = (3) 
In Fourier space, the power-spectrum of the noise is given by 

/+oo /<+oo 
dtexp(iojt)S(t) = dtexp(iwt)(£(t)£(0)}. (4) 
-oo J — oo 

If £(t) is a white noise of amplitude T>, we have 

S(t' -t) = VS(t' - 1) and S{w) = V. (5) 
When is an Ornstein-Uhlenbeck process of amplitude T> and of auto-correlation time r, we have 

S(f-t) = £e-l*- t 'l/ T and = — ^ . (6) 

2 r 1 + w'r 

In this work, we shall consider the case where the power-spectrum of decays at high frequencies in the following 
manner 

S(u) ~ V{ujt)- 2,j when |w|->oo. (7) 

The amplitude T> of the noise and the correlation-time r arc defined by dimensional analogy with equation © . The 
exponent a characterizes the high frequency behaviour of the power-spectrum. When a is an integer, such a noise 
can be generated from the white noise by solving a linear differential equation of order a . 

We shall prove that in the long time limit, the statistical properties of the oscillator in the phase space can be 
classified by the following two parameters (i) the exponent v, defined in Eq. that encodes the large amplitude 
behaviour of the confining potential U\ (ii) the exponent a that determines the high frequency behaviour of the 
power-spectrum. For fixed values of v and a, the phase space distribution takes in the long time limit a universal 
form (that also depends on the dimensional parameters 7, r and T>) that we shall calculate. 



II. EFFECTIVE DYNAMICS IN THE ASYMPTOTIC REGIME 



A. Use of integrability 

The mechanical energy of the oscillator is defined as 

E =^x 2 +U(x). (8) 

In the absence of noise and dissipation this quantity is conserved. This implies that the Hamiltonian system underlying 
Eq. ([T]) is integrable. It is therefore possible to define an action variable J and an angular variable (f> so that the 
transformation (j> = x, x) — > (J, 4>) is a canonical transformation. For a given value E of the energy, the angle 4> is 
given by !EH 

<j) = u(e) r . dy w^ oj{e) = (~ r max - dy -\ , (9) 



y/2(E-U{y)) y/2{E -U{y)) 

where x max satisfies U(x max ) = E (recall that Li is an even function). With this definition, the range of the phase <j) 
is 2ir. The action variable J is a function of the energy only and is determined by the following equation: 

The Hamiltonian equations of motions, in terms of the action-angle variables, read simply: 

j=-^iM = (11) 
dep 

■ dE(J,<f>) 

cj) = — — — =w(E). (12) 



The variables (E, <p) define a bona fide set of coordinates on phase space. The formulae for transforming the variables 
from position and velocity to energy and angle are given by: 



+00 



x(E,4>) = J2 x n{E)j n * (13) 

n— — 00 

+00 +00 

±{E,4>) = J2 ME)e m * = ME) J2 ™ n (E)e in *. (14) 

n— — 00 n— — 00 

Here again to write this change of variables, we have used only the deterministic and dissipationless parts of the 
dynamics. We choose the origin of 4> such that x n (E) is a real number and that x n (E) = X- n (E). We also write 

\x 2 {E,J>) = £ y n (E)e in +. (15) 

71 — — OO 

We then have 

x / n— — 00 

- \-x\EA)) = ME) E nVn(E)e in +. (17) 



n— — 00 



In the general case, we find from Eq. JT]) that the time variation of the energy is given by 

dF = -^ (18) 

We have used here the rules of classical calculus when changing variables 0, H^|: this is allowed because we are 
working with the Stratonovich interpretation of Eq. ([1]). The energy variation has thus two contributions: a loss term 
due to friction and a stochastic 'elastic energy' term due to the work of the random multiplicative force x£ (t) . 
The time variation of the angle variable is given by 

d(j) , r z. ~.rfu.\\ fn\® x fn<\ (i?\z.® x ^\ 3 (1 2 



]t = U j(E)M.ji-xm)u(E)^=uj(E)+ 1 u(E)x—-atME)—l-x 2 (E, ( f>)) . (19) 



If we substitute in the two dynamical equations (|18|) and (|19|) . the expressions given in Eqs. (|13|) . (|14j) . (I16|) and 
(fTT)) . we obtain 



+00 +00 +00 



£ = 7 w 2 (£) E z»(£)z m (£)e i(n+m)0 -£(iM£) E iny n (E)j n * (20) 

n=— 00 m— — 00 n— — 00 

= W (£0+V(^0 £ E in*n(£0^^e*(^ + e(t)a;(£0 £ %^e^ . (21) 



n— — 00 m— — 00 n— — 00 



We emphasize that this coupled system of stochastic nonlinear equations is rigorously equivalent to the initial random 
dynamical equation ([1]). 



B. The effective Markovian description 

Although the problem we study here is non-Markovian, because the noise has a non- vanishing correlation time, it is 
possible to derive for the associated Probability Distribution Function Pt(E, <p) an effective coarse-grained Markovian 
equation using a Kramers-Moyal type expansion @, H, [26[ : 

n— 1 m+fc=n x 7 x 7 

m,fc>0 



where we have defined 

M m}k (E,4,,t,S) = ((AE t (S)) m (AM5)) k ) , (23) 

with AE t {8) = E(t + 8)-E(t) = I dsE(E(t + a),4>(t + s),t + s) , (24) 

Jo 

and A(f> t (5) = cj>(t + 5) - 4>{t) = [ ds<j>(E{t + s),(j)(t + s),t + s) . (25) 



The expressions of E and 4> on the r.h.s. of Eqs. f2"4"|) and (|25p are given in Eqs (|20p and (j2"Tj) respectively. 

The time scale 6 that appears in the Kramers-Moyal expansion must be chosen in a physically relevant manner: 
8 has to be small but must remain larger than the intrinsic period of the oscillator (this condition is automatically 
fulfilled at large amplitudes because the intrinsic period tends to zero). Besides, one must also have S ^S> r (where 
r characterizes the correlation-time of the noise) in order to end up with an effective Markovian description of the 
dynamics. 

A systematic procedure for evaluating the coefficients M m .fc that appear in Kramers-Moyal expansion has been 
developed by Carmeli and Nitzan in [24]]. We first rewrite Eqs. (j2"4f and ([2"5"|) as 

AE t {5) = [ dsE(E(t) + AE t (s),4>(t)+A4> t (s),t + s) (26) 
Jo 

A&(J) = / ds <j> (E(t)+ AE t (s), <f>(t) + A<f> t (s), t + s) . (27) 
Jo 

The values of AE t (6) and A<fi t (S) are evaluated according to the following iteration scheme, labelled by the integer 
index I 

AEf' (8) = J dsE(E(t)+AE( l - 1 \s),<p(t)+A4 l ~ 1) (s),t + s) (28) 
A4 l - 1] (6) = ds^E^ + AEt'Hs^^+A^is^t + s) , (29) 
with initial values given by 

A£ t (0) (s)=0 and A^ 0) (s) = u>(E)s . (30) 

Performing this expansion and neglecting the terms of order S n with n > 1, we obtain, in the limit S — > 0, after 
some systematic but tedious calculations: 

\ { AE { 8)) = -W(E) £ n>\ Xn (E)f + £ ^ ( + ^ {E) ^gll!} 

+ ^ E - 2 \yn(m 2 %, (3D 

n— — oo 

j- § ((AE)*(6)) = ^ g n 2 | y „(i?)| 2 5„, (32) 

n— — oo 

(AE(6)A<I>(5)) = 0, (33) 
where we have defined 

/+oo 
dtexp(inw(E)t){£(t)£{0)) . (34) 
-CJO 

We do not need to give the exact values of (A0(<5)) and ((Acj)) 2 (5)) because they will have no incidence in the following 
calculations. Higher moments are negligible at the considered order of the calculations. 



C. Effective Fokker-Planck equation for the energy 



We now substitute the average values calculated in Eqs. (|31I32I and [3"3"|) into the Kramers-Moyal expansion 
Because the cross-correlation term (AE(6)A<f>(6)) vanishes, we can integrate out the angular variable from Eq. 
and obtain an effective Fokker-Planck equation for the energy: 



dt dE { 6 " K ' J dE 2 1 25 

Defining the following two auxiliary functions: 



+00 



ex(£0 = J2 n 2 \x n {E)\\ (36) 



+ OO 



^(E) = 5 £ n 2 | y „(i?)| 2 5„, (37) 



2 

n— — 00 

we rewrite Eqs. (pH"]) and as follows: 

~<A£(«5)) = - 1 u; 2 (E)e 1 (E)+u; 2 (E)^l+e 2 (EME)^l : (38) 
1((AS) 2 (<5)) = ^(i?)^^). (39) 
Substituting these expressions in Eq. (|3"5)) leads us to the effective Fokker-Planck equation for the energy: 

^p- = ^^(E)(je 1 (E)+e 2 (E)-^jcj(E)P t (E)^ . (40) 

For dissipationless motion, 7 = 0, this equation does not have a stationary solution: the particle diffuses in phase 
space by absorbing energy from the noise and there is no mechanism to limit the growth of the amplitude. The 
observables grow as power-laws with time as the explicit solutions of the next section will show. When 7 7^ 0, the 
system reaches a stationary measure characterized by a non-equilibrium steady state with an asymptotic probability 
distribution fgtatC^) that differs from the canonical Boltzmann-Gibbs law. For 71 < 1, the effective Markovian 
description remains valid and the stationary solution of the effective Fokker-Planck equation is given by 

W^ = ^e^ d «« (41) 



where the prefactor J\f ensures the normalization of P s tat C^) ■ 



III. EXPLICIT SOLUTIONS 



A. The Hamiltonian case 



We shall first consider the case where the dissipation effects are not taken into account. In the absence of dissipation, 
the physical observables such as the amplitude, the velocity and the energy of the oscillator grow as power-laws with 
time. We shall calculate the associated scaling exponents and prove that their values depend only on v which 
determines the behaviour of the external potential at large amplitudes and on a that measures the relative weight of 
high frequencies in the noise spectrum (and which also characterizes the smoothness of the noise) . 

In the long time limit, the particle diffuses to large amplitudes in phase space. Therefore in Eq. {T]) we can neglect 
the linear restoring force (proportional to luq) and replace the potential by its asymptotic behaviour given in Eq. 

U{x) ~ ^— . Then, the change of variables to energy and angle coordinates, given in Eqs. (JSJ and ©, takes the 
simpler form: 



The equation of motion for the underlying deterministic system are given by 

E = (43) 

4> = u(E) where u(E) = / - = C V E^ with C„ = ^ " f^ V^, (44) 

where the last formula, in terms of the Euler Gamma function, is obtained from |40]. We now define [4l[ the 
hyperelliptic function T v : 



' ' (45) 



%(Y)=X «->F = / 
Jo 

The function T v is periodic with period 

rl du _ 2^ r(^) 

/ Vi - ^ «> r(^) 

Inverting Eq. (|42p . we express the position and the velocity in terms of energy and angle using the function T v pol. Ell] : 



y/1 - n 2 " 



^=4/ ^^^^T- (46) 



x(£,0) = (2vE)^T v (^ , (47) 

£(25,0 = V2ETl(^j , (48) 

where 7^ is the derivative of the function % which, using Eq. (|4*5|) , satisfies the relation 

(r,(y))^ + (r;(r)) 2 = i. (49) 

The coordinates x and x are 27r periodic functions of the angle variable (j>; they can thus be developed into Fourier 
Series as in Eqs. (fl3|) and (fl4l . More precisely, if we write 

^(r^)= £ /nem0 ' (5o) 



n— — oo 



we obtain the Fourier coefficients of x and x 

x n (E) = [2vE)^f n and v n {E) = V2E (m/„) . (51) 

We note that, in the present case, the Fourier coefficients depend on the energy E only through a global prefactor 
that does not depend on the harmonic n. 

This identification allows to calculate exactly the function e±(E) defined in Eq. 



n— — oo 



2 



= ^4^\ W ^^ dy= (H^jf 1 d „ V T^, (52, 

where the second equality is obtained using Parseval's identity, the third equality using Eq. (|49[) over a quarter of 

a period of the hyperellectic function T v and the fourth equality by the change of variable T V {Y) = u. Using the 
expression (|46|) and evaluating the last integral in terms of Gamma functions [401 ] , we obtain: 

, , (2uE)iv T 2 (^) , 

We now calculate the function 62(E), defined in Eq. (|37[) . Using Eqs. (fT"5)) and (ITT)) , we find that 

yn (E) = { ^p^g n where = g 5 „e^ . (54) 



Because T v is a real and even function of (f>, we have g n = g- n and we can rewrite 62(E) as follows: 

e2 (E)= { -^Y. n2 9nS n . (55) 

71=1 

This sum depends in a non-trivial manner on the noise spectrum. We are, however, interested in the long time 
behaviour of the Probability Distribution Function dPt (E) . When t — > 00 the typical value of the energy E of the 
system also increases without bounds and, therefore, the intrinsic frequency tu(E) of the system, which according to 
Eq. (TS)) is proportional to E^~ , also increases without bounds. Thus, when t — > 00, we can replace <S„( = S(nu(E)) ) 
by its asymptotic behaviour given in Eq. (J7J and obtain 

x v 7 J n—1 

Denoting by A a the value of the convergent series J2n=i n2 ~ 2 ° 9n? we can write 

2 

t2(E) = '^E 2 -^ 1 W i th Vv<a=V WL. Aai (57 ) 

where C„ was defined in Eq. (fH]) . 

When the noise is white a = 0, its power spectrum is constant and the sum in Eq. (|55|) can be evaluated exactly. 
Following steps similar to those which led to Eq. ([5^]) . we obtain 



s (2vE)% r(£)r(&) 
w r(^)r(^ 



For an Ornstein-Uhlenbeck noise a — 1 and we must evaluate the expression Xm=iSn' This, again can be done 
explicitely, thanks to the Parseval identity: 



, Mg f r(^)r(A) r^) 
w Vr(^)r(^) 1*0*3; 



We now deduce the asymptotic expression of the Probability Distribution Function dPt (E) in the limit t — * 00 and 
when there is no dissipation. The effective Fokker-Planck equation then reduces to 

dPt{E) - 9 LiEME) 8 "™™) ■ (60) 



dt as V dE 

Taking into account the expressions of uj(E) and 62(E), given in Eqs. ([331) and ([57)) respectively, we can rewrite this 
equation as 

^^AC,?d( E ,dE^P t (E)\ ^^^ + 3-2^-1) 

dt t 2 ° dE y dE J 2v v ' 

This equation has a self-similar structure [3^ | and it is natural to look for solutions of the form 

*(*)-i*(-fl) w,«, and a-fe±^, (62) 

the prefactor 1/E ensures that Pt(E) is normalized. The function <j>(u) of the scaling variable u = E a /(Kt) satisfies 
an ordinary differential equation. The solution of this equation is given by 

(j>(u) oc u^e~^ . (63) 

Inserting this solution into the expression (]61|) for P t (E) we find, after normalisation, the following asymptotic formula 
for the probability distribution function: 

a 1 / E a \^ ( E a \ 



where K and a are defined in Eq. (|62p . 

From this general result, we can retrieve the solutions for white noise (W) and for Ornstein-Uhlenbeck (OU) noise. 
For white noise, using Eq. (|58|) we have 



1 f E z=l}*F=V C E i 

,•/*-•> — . 



P™(E) 



r(-f+i ) vE \2V w tl 1 22? w t 



Wlth Pw " P 2^7+1) r ( JL)r ( ^)- (65) 

This formula was also obtained in [l7| by stochastic averaging. For Ornstein-Uhlenbeck noise, we have, using Eq 

. "+ 1 



ou, ^ 1 2(2/ -1) E^\^ E^ 



r 



v+i \ vE 

4(^-1), 



with v v f"- 1 !'^ ^ 1 £!ffiE!iil^ (66) 

This formula is identical to the one derived in [l8| by using the averaging method at the second order. 

To summarize, we have shown that in the long time limit the scaling behaviour of the dissipationless nonlinear 
oscillator in presence of a noise with a power-spectrum that satisfies equation j7]), is given by 



E 



x 



Vt \ 



Vt \ 2 ("+ 1 )("- 1 ) 



/ T)f\ 2(<r + l)(„-l) 

These scalings, derived here by a systematic calculation, agree with the results given in 20] that were conjectured 
by performing a partial resummation of the small correlation-time expansion of the stochastic Liouville equation 
associated with the random dynamical system under study. The method used in (T5 , H3] was an approximation that 
could not be applied to the system studied here but only to a simplified model. In fact the resummation technique 
yielded the correct exponents but the prefactors were out of reach. Here, we have obtained the closed expression 
for the Probability Distribution Function, which contains the full information on the statistics of the system when 
t — ► oo, i.e., it provides us the scaling exponents and the corresponding prefactors. 



B. The dissipative case 



In the presence of dissipation, the system reaches in the long time limit a steady state with the stationary probability 
given by Eq. pTj) . Using Eqs. ([55]) and (T55|) . we find the explicit formula for this stationary distribution: 



a 1 / 2vy E a \ *™ ( 2vj E a 



p "W = ^Ei{ricw) ex H-^Ti)J' (68) 

where K and a are defined in Eq. (j6"2"|) . From this stationary PDF we find that the typical value of the energy is given 
by 

Ea ~ ^ ■ (69) 

We observe that the expression of P s t a t (E) becomes identical to that given in Eq. ([6^)1 if the time t in Eq. ([6^)1 is 
replaced by t 7 with 



The value of t 7 determines the time scale at which the system becomes sensitive to the dissipative effects. For 
(Ci 7 , the system evolves as if it were Hamiltonian and the physical observables grow algebraically with time. For 
t > i 7 , the system sets in its steady state and the statistical averages of the physical observables become stationary. 
Finally, we justify the validity range of the stationary probability distribution ([55)1 . We recall [2J] that the Markovian 
approximation is valid only if 7 -1 is much smaller than the correlation time r of the noise, i.e., 

7T<1. (71) 

Besides, the presence of dissipation should not alter significatively the dynamics of the fast variable 4>: thus on the 
r.h.s. of Eq. (f2Tj) the second term must remain much smaller than the first one. Using Eqs. (|4"4"1 H71 and I38]l this 
requires that 7 <C E^^ 1 ^ 2u . From the typical value (|69ll of the energy and the expression of a given in (|62ll . this 
condition becomes 

(7t) 2ct 7 3 < V ■ (72) 

Thus, a sufficient condition is 7 s < T>. For higher values of the dissipation rate, the expression will no more be 
valid. The behaviour of the system can change drastically and a phase transition to a state localized at the origin 
x — x — can occur 0, [HI, [H| ■ 



C. The additive noise case 



In this last subsection, we study the nonlinear oscillator driven by an additive noise: 

dF»(*)+7^(t) + -^=«t). (73) 

The change of coordinates to energy and angle variables is the same as in Eqs. (fTS"]) and (JHJ). The dynamical equations 
in these coordinates will read 

+00 +00 +00 
E = 7^ 2 (£) J2 J2 Xn{E)x m (E)S n+m ^ + £(t)w(E) inx n (E)e in * (74) 

n— — oQ m= — 00 7t— — 00 

n— — 00 m— — 00 n— — 00 

Using the Carmeli-Nitzan technique, we derive the following effective Fokker-Planck equation for the energy: 

^^ = ^{u{E)(^e 1 {E)+Um^u{E)P t {E)} , (76) 
the function t\(E) was defined in Eq. (j3"6"|) and 62(E) is given by 

UE) = \ £ n 2 |x„(£;)| 2 5„. (77) 

n— — 00 

Denoting by A a the value of the convergent series Yln=i n2_2fT /n> we can wr ite 



'AE^^E^ 1 with Vur-V-^X. CM 



The constant C v was defined in Eq. (|44|) . It is possible to carry out explicit calculations following the same lines as 
for the multiplicative noise case. If dissipation is neglected, probability distribution function is found to be 



a 1 f E a \ ^ f E° 



where a and K are given by 



«< B > = vm e K&Tt) "l-SiR.I- (79) 



, = ( a + 1) („- 1) + 1 md R = v^^ 



For the special cases of white noise or Ornstein-Uhlenbeck noise, explicit expressions for K can be derived and the 
formulae for the probability distribution function derived in [l9l ] using the averaging method are recovered. 

Thus, in the absence of dissipation, the following algebraic scalings for the main observables of the system are 
satisfied 



E 



Dl \ („+i)(„-i)+i 

_2ct 



Dl ^ 2(,x + l)(„-l) + 2 
_2cr 



T>t \ 2(<r + l)(^-l) + 2 



r'2a 



(81) 



In particular, when the noise is white (or if the potential is quadratic), we recover the result that the energy grows 
linearly with time. The amplitude of the noise being constant, the exponents in the additive case are smaller than 
those in the multiplicative case, as expected. 

If we take dissipation into account, the system reaches a steady state in the long time limit. The expression of the 
stationary probability matches that of Eq. (fT5|) if t is taken to be 

^72(^ + 4-1))' (82) 
This expression defines the dissipation time scale for a nonlinear oscillator subject to additive noise. 



IV. CONCLUSION 



We have used an effective coarse-grained Markovian description to carry out a quantitative analysis of the nonlinear 
oscillator confined by a polynomial potential and subject to Gaussian noise of arbitrary spectrum that decays as a 
power-law at high frequencies. This approach has allowed us to calculate the distribution, in phase space, of the 
dynamical system in the long time limit. In the absence of dissipation, the particle diffuses without bounds with 
an anomalous scaling law. The energy transfer from the random perturbation to the particle also follows a scaling 
power-law. The diffusion exponents and the corresponding amplitudes are determined exactly. In the presence of 
dissipation, the system reaches a non-equilibrium steady state; the corresponding stationary distribution has also 
been calculated analytically in the limit of vanishingly small dissipation. The advantage of the method used here as 
compared to stochastic averaging is that it can readily be adapted to any noise spectrum. For white and Ornstein- 
Uhlenbeck noises, the two approaches give identical results. It would be of great interest to apply this method to 
higher dimensional integrable systems subject to stochastic perturbations such as the nonlinear Schrodinger equation 
in a random potential, which is often used to study the effect of nonlinearity on localization. 
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